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ABSTRACT 

We present the detection of infall, rotation and outflow kinematic signatures towards both a protostellar 
source, VLA 1623 and what was initially thought to be a pre-protostellar core, SM1N, in the p Ophiuchus A 
region. The kinematic signatures of early star formation were detected in the dense molecular gas surrounding 
the embedded sources using high signal-to-noise millimeter and submillimeter data. Centroid velocity maps 
made with HCO + J=4— >3 and J= 1 — >0 line emission exhibit the blue bulge signature of infall, which is predicted 
to be seen when infall motion dominates over rotational motion. Further evidence for infalling gas is found in 
the HCO + blue asymmetric line profiles and red asymmetric opacity profiles. We also performed CO J=3— >2 
and J=l — >0 observations to determine the direction, orientation, and extent of molecular outflows, and report 
the discovery of a new bipolar outflow possibly driven by SM1N. 
Subject headings: stars: circumstellar matter-ISM: clouds-stars: formation 



1. INTRODUCTION 

Observationally, the earliest stage of evolution for low-mass 
protostars is the so-called "Class 0" phase, which is character- 
ized by a blackbody-like spect ral energy di stribution (SED) 
that peaks in the submillimeter (Andre et al. ( 1993), hereafter 
AWB93). Barsony (1994) defines the following as charac- 
teristics of Class objects: Lboi/L\^ mm < 2 x 10 4 , an SED 
like a 30 K blackbody, undetected at A < 10/im, and posses- 
sion of a molecular outflow. These characteristics of Class 
objects imply that they are in a very active mass accre- 
tion phase, and possess more circumstellar mass than stellar 
mass. An even earlier phase of evolution prior to class ob- 
jects h as also been proposed, the so -called "pre-protostellar" 
stage (Ward-Thompson et al. 1994). Observationally, objects 
in this stage of evolution have SEDs that are similar to Class 
objects, but they lack other kinematic signatures of star 
formation such as molecular outflows and centimeter con- 
tinuum emission. The pre-protostellar cores are probably 
gravitationally bound, but lack a central hydrostatic stellar 
source, and are hence not yet "protostars". Indeed molecular 
line spectroscopy o f a number of pre-protostellar cores (e.g. 
Tafal laet alJ fT998)) show that they are undergoing infall, but 
more slowly than that seen towards Class protostars. 

To study the detailed kinematical evolution of pre- 
protostellar objects to Class objects and beyond, it is impor- 
tant to study both type of objects with observational tools like 
millimeter and submillimeter spectroscopy. Pre-protostellar 
cores do not possess a strong outflow, and hence it is ex- 
pected to be less complex to disentangle the various motions 
that accompany star formation such as infall, rotation, out- 
flow and turbulence. Since the details of the star formation 
process can vary from cloud to cloud, and may depend on 
environmental factors, it is crucial to choose pre-protostellar 
and protostellar candidates within the same region. The cen- 
trally condensed L1688 core of the nearby (160 pc) p Ophi- 
uchus (hereafter p Oph) molecular cloud consists of a rich 
array of young stellar o bjects at various stages of evolution 
( Wilking & Ladalfl983l) . In particular, the cloud A region 
(called p Oph A) harbors a number of pre-protostellar and 
protostellar cores. In this region, there is VLA 1623, which 
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was designated as the prototype of the Class phase of evo- 
lution (AWB93). VLA 1623 is known to possess a remark- 
ably colli mated and very young (~ 6000 yrs) molecular out- 
flow (lAndre et all (rT9"90). AWB93). The millimeter and sub- 
millimeter continuum maps of AWB93 reveal that p Oph A 
breaks up into 4 cores, SMI, VLA 1623, SM1N, and SM2. Of 
these, VLA 1623 was the only source with an identified out- 
flow. Both SM1N and SMI are probably real pre-protostellar 
cores (A WB93). Wide-field 1.3 mm continuum mapping by 
(Mott e et all 19981) (hereafter MAN98) reveals that Oph A is 
the brightest of all Oph cores, and in addition to VLA 1623 
there are eight other starless cores in p Oph A. The objects A- 
MM6, A-MM7, and A-MM8, together with SM1N, SMI and 
SM2, compose an arc within Oph A, which seems to be de- 
lineating a photo-disassociation front of a PDR region illumi- 
nated by the nearby B star SI (see MAN98). Recent high an- 
gular resolution N2H" 1 " observations (Di Francesco et al. 2004) 
of the quiescent dense gas in this region show emission peaks 
of this tracer at SMI and SM1N, but not at VLA 1623. N 2 H + 
can deplete signifi cantly at densities greater than 10 5 cm" 3 
(Bergin et al. 2002), so the difference in N2H + emission from 
SM1N to VLA 1623 might be another indicator that SM1N is 
less evolved than VLA 1623. 

In this paper, we present results of a study of the p Oph 
A region designed to make a comparison of pre-protostellar 
core(s) with a more evolved Class protostellar object 
(VLA 1623). One of the main goals of this study was to under- 
stand better the earliest stages of star formation from the pre- 
protostellar phase to the Class phase. However, as will be 
seen from the results of this study, SM1N seems to power an 
outflow, and is hence probably not pre-protostellar. To probe 
the kinematics and physical conditions of the cores, we per- 
formed a multi-transitional study using both density sensitive 
molecular tracers, and CO and its isotopes. HCO + and N2H + 
preferentially trace high density gas, while CO traces a much 
wider range of gas densities, in addition to being a good tracer 
of outflows. Submillimeter transitions of HCO + and CO probe 
regions of higher density and excitation than the millimeter 
transitions. For this reason, we augment our millimeter ob- 
servations with submillimeter observations of corresponding 
transitions. In §2, we describe our observations, and in §3 we 
present our results. 
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TABLE 1 
Observations 



Molecule Transition Observatory Frequency (GHz) Beam Width (") 



H 13 CO + 


J= 1 


-0 


FCRAO 


86.754288 


62 


HCO+ 


J= 1 


->0 


FCRAO 


89.188496 


60 


N^H + 


123^ 


012 


FCRAO 


93.173770 


58 


13 CO 


J= 1 


->0 


FCRAO 


110.201354 


49 


CO 


J= 1 


-+0 


FCRAO 


115.271202 


47 


13 co 


J = 3 


->2 


CSO 


330.587980 


23 


c 17 o 


J = 3 


->2 


CSO 


337.061123 


22 


CO 


J = 3 


->2 


CSO 


345.795990 


22 


H 13 CO + 


J = 4 


^3 


CSO 


346.998338 


21 


HCO+ 


1 = 4 


^3 


CSO 


356.734134 


21 



2. OBSERVATIONS 

2.1. FCRAO 

Millimeter observations of HCO + , H 13 CO + and N 2 H + to- 
wards VLA 1623 and SM1N were performed with the Five 
College Radio Astronomy Observatory (FCRAO) 14m tele- 
scope in December 1998, February 1999 , and June 2001 using 
the SEQUOIA 16-beam array receiver (Erickson et al. 1999), 
and the FAAS backend of 16 autocorrelation spectrometers. 
Mapping of the 12 CO and 13 CO J=l >0 transitions were also 
performed in February 2004 using the newly-developed on- 
the-fly (OTF) mapping technique at FCRAO using the up- 
graded SEQUOIA array receiver which had dual-polarization 
capability, and a dual channel correlator (DCC) backend that 
allowed the observation of both the 12 CO and 13 CO isotopes 
simultaneously with a spectral resolution of 50 kHz. The 
2004 observations covered a 12' x 12' area of the p Oph 
A cloud. Table ^ summarizes the frequencies and effective 
resolutions of each transition. HCO + and N2H" 1 " line transi- 
tions were observed via frequency-switching with subsequent 
folding and third-order baseline removal. HCO + and its iso- 
topic observations span a region 6' x 6' centered on VLA 1623 
(a(1950) = 16 /, 23" , 24.8 I ,<5(1950) = -24°17'46"). Pre-OTF 
data were reduced and analyzed with the Gildas CLASS soft- 
ware package. The OTF data was reduced using the FCRAO 
OTFTOOL software package (iHever etalj r2004). Pointing 
and focus were verified with observations of SiO masers. 

2.2. CSO 

We obtained sub-millimeter observations using the 10.4 m 
telescope of the Caltech Submillimeter Observatory (CSO) at 
Mauna Kea, Hawaii. The May 1996 HCO + and CO obser- 
vations were performed using On-The-Fly (OTF) mapping. 
These observati ons were made wit h the 345 GHz SIS waveg- 
uide receiver JEllison etal.i n"989) and a 1000 channel, 50 
MHz wide acousto-optical spectrometer. HCO + observations 
span a region 70" x 70" centered on both VLA 1623 and 
SM1N. Our CO observations mapped a region of 150" x 150" 
centered on VLA 1623 and SM1N. Individual OTF maps 
made on VLA 1623 and SM1N were combined together to 
derive the results presented below. Data were reduced and 
analyzed using the CLASS software package. 

3. RESULTS 

3.1. Molecular Line Profiles 

Figure ^ shows a summary of the spectral observations 
towards the central positions of two sources, VLA 1623 
and SM1N. The left panel shows the observations towards 



VLA 1623 and the right panel shows the observations towards 
the source SM1N. Table|5]lists the derived centroid velocities 
for various molecular transitions towards the central position 
of VLA 1623 and SM1N. In order to exclude any effects of 
outflows, the centroid velocities in Table |2] were computed 
within a velocity interval of 2 to 6 km s" 1 . 



(a) (b) 




FIG. 1. — Observed line profiles toward the central positions of (a) 
VLA 1623 (left) and (b) SM1N (right). The multiplication factor is stated 
after the names of each line. Both velocity ranges are -5 to 13 kms~'. The 
vlsr of ~ 3.6 kms~' obtained using N2pf + is shown by the vertical dotted 
lines. The systemic velocity of the cloud is ~ 4 km s . 

Since N2H + is typicall y optically thin, and is a good tracer 
of dense, quiescent gas (Benson et al. ( 1998) and references 
therein), we used the J=l — >0 observations of N2H + to deter- 
mine the true systemic velocity tow ards these sources . Hyper- 
fine fitting of the N 2 H + lines (c.f. lCaselli et alJI 19951) yields a 
\lsr of 3.63 ± 0.02 and 3.58 ± 0.01 km s _1 towards VLA 1623 
and SM1N respectively. 

The central spectra of HCO + and CO towards VLA 1623 
and SM1N show strong self-reversed dips. These dips imply 
true self-absorption as the optically thin tracers (N 2 H + , C 17 0, 
H 13 CO + ) peak at the velocity of the self-reversal. One inter- 
esting feature of the submillimeter CSO HCO + J=4^3 spec- 
tra is the lack of the "classic" blue asymmetry expected in in- 
falling regions. It is this "classic" blue asymmetric line profile 
signature that has been the technique of choice f or mo st stud- 
ies of infall (see for example, Gregers en et alJ ([11973, here- 
after GEZC97). Indeed, the lack of blue asymmetry in their 
HCO + observations led GEZC97 to conclude that VLA 1623 
was an ambiguous case for the interpretation of infall. In fact, 
it can be seen from Figure ^that the central HCO + J=4^3 
line profiles in both sources exhibit a red asymmetry, a pro- 
file consistent with expanding rather than infalling gas. In- 
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TABLE 2 

Centroid Velocities of Central Line 
Profiles 



Source 


Transition 




v c (kms ') 


vlsr (kms -1 ) 








VLA 1623 


CO J=l— >0 




3.70 ±0.02 


(3.6 a ) 


l3 CO J=l — >0 




3.36±0.01 




CO J=3— >2 




4.01 ± 0.01 




13 CO J=3^2 




3.32 ±0.03 




C 17 J=3^2 




3.53 ±0.06 




HCO+ J=l-»0 




3.06 ±0.36 




H CO J=l — i 





3.yi ±0.11 




HCO+ J=4^3 




3.56 ±0.03 




H 13 CO+ J=4— 


3 


3.81 ±0.06 


SM1N 


CO J=1^0 




3.60 ±0.02 


(3.6 a ) 


13 COJ=1^0 




3.32±0.01 




CO J=3^2 




3.60±0.01 




13 CO J=3^2 




3. 16 ±0.04 




HCO + J=l-»0 




2.9 ±0.4 




H 13 CO + J=I — 





3.34±0.16 




HCO+ J=4^3 




3.53 ± 0.04 




H 13 CO + J=4— 


3 


3.75 ±0.09 



Via derived from NjH* hyperfine structure 



deed, a very strong and collimated outflow has been discov- 
ered towards VLA 1623 JAndre et alj n"990). It is likely that 
the expected blue-asymmetric line profile signature of infall 
has been inverted in the presence of this strong outflow. The 
isotopic H I3 CO + J=4^3 has a centroid velocity (see Tabled) 
which is red-shifted with respect to the main line (as might 
be expected for infall), but at the same time is slightly more 
red-shifted with respect to the ~ 3.6km s -1 velocity derived 
from the N2H + transition. In the presence of pure infall, the 
optically thin H I3 CO + line would be expected to be centered 
on the sy stemic velocity of the cloud core. As will be shown 
in 33.2.21 HCO + and its isotopes do trace the outflows, so this 
may cause the velocity centroid of these lines to move sys- 
tematically redward, compared to the pure infall scenario. 

In comparison with their corresponding HCO + J=4^3 
transitions, the HCO + J=l— >0 spectra are heavily blue 
shifted. Their self-absorption dips almost reach the contin- 
uum baseline, and the red-shifted peak is barely visible within 
the signal-to-noise ratio of our observations. Why are the 
HCO + J=l — >0 observations so different in appearance from 
the HCO + J=4^3 observations? The two transitions are 
likely probing differing regions of gas. The gas observed by 
the larger beam of the FCRAO 14 m in HCO + J= 1 — >0 comes 
from extended regions around the two sources. The extended 
envelope traced by the lower transition may be indicative of 
the overall large scale motion of the gas due to collapse. The 
CSO HCO + J=4^3 transition probes a smaller region of gas 
as well as probing farther into the envelope. Thus, the shorter 
wavelength submillimeter transition may be picking out the 
densest gas regions near the protostellar cores. Probing nearer 
the source allows detection of the inner region of infall but 
also the developing outflow, which might cause the red asym- 
metry seen in the submillimeter transitions. 

The dipole moment of CO is much smaller than HCO + and 
is thus expected to trace a wider range of densities. CO is 
often seen tracing swept-up and outflowing material, whereas 
the HCO + isotopes probe more of the denser, centrally con- 
densed gas. Indeed the broad wings of emission seen for 
the CO J=3^2 spectra in Figure [2 imply that CO is tracing 



outflows in both sources. Curiously, the central line profiles 
of CO also show the "classic" blue asymmetry already dis- 
cussed. The CO spectra show deep self absorption in the line 
core and very broad wings of emission that extend beyond 
±10 kms -1 from the systemic velocity. The self-absorption 
dip and blue asymmetry profile of the I3 CO J=3— >2 support 
the notion that this isotope is optically thick as well. The rarer 
isotope C 17 has a gaussian line profile which peaks at the 
Vlsr of VLA 1623, as expected for an optically thin line. 

In the presence of infall, the line emission from the rarer 
isotope is expecte d to be r ed-shifted with respect to that from 
the main isotope (Naravanan et al. 19981 12002I) . In the cases 
of VLA 1623 and SM1N, all of the centroid velocities of the 
rarer isotopic spectra are red-shifted with respect to their main 
lines (see Table [2j. If infall were the only source of mo- 
tion, both the emergent line intensity and optical depth profile 
must be asymmetric, with the line profile showing blue asym- 
metry, and the optical d epth profile showing red asymmetry 
(Nara vanan et aljri998l) . In the presence of pure expansion 
(outflow), the converse would be true: the line profile would 
be red asymmetric, and the optical depth profile would be blue 
asymmetric. 

We performed an optical depth analysis using the CSO 
HCO + and H 13 CO + J=4^3 observations. In Figure we 
present the optical depth profiles for VLA 1623 and SM1N. 
The optical depth calculations were performed only over the 
FWHM of the isotopic lines. The opacities were estimated 
from the observed line profile ratios of HCO + /H I3 CO + . The 
opacity, r^f of the more abundant isotope was estimated us- 
ing I(Main)/I(Iso)= (1 -e T *)/(l -<? r * '), where r is the ratio 
of main to isotopic abundances, and I(Main) and I(Iso) are 
the line intensities of the main and isotopic lines respectively. 
An isotopic abundance ratio of [HCO + ]:[H 13 CO + ]= 45:1 was 
used. This technique assumes that the excitation temperature 
is the same for the main and isotopic species. 



(a) (b) 




FIG. 2. — HCO + J=4— >3 optical depth profiles toward the center positions 
of (a) VLA1623 (left) and (b) SM1N (right). The temperature scale is shown 
on the left-hand side of the Y axis, while the opacity scale shown on the right- 
hand side Y axis. The velocity for the ambient cloud (3.6 kms -1 ) derived 
from our N2H + observations is shown as the vertical lines. The line profile of 
HCO + is shown in the solid red histograms and the opacity profile is shown as 
the dashed blue histograms and only over the velocity range of the H 13 CO+ 
line profile. 

From Figure |2 we see that not only are the main line pro- 
files red asymmetric towards both VLA 1623 and SM1N, but 
the optical depth profiles also exhibit a red asymmetry. The 
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red-shifted line profiles are consistent with expansion or out- 
flow but the red-shifted opacity profiles are consistent with 
tracing infall. The main isotope HCO + J=4— >3 must be trac- 
ing a significant part of the outflow for its line profile to show 
an expansion signature. However, the opacity profile which 
is principally determined by the more optically thin isotope 
seems to be more consistent with tracing infall. The optically 
thin isotope is likely tracing material further towards the cen- 
ter of the cloud cores, whereas the main line is likely being 
contaminated by the expanding (dense) shells of the outflow. 
This indicates the importance of using both isotopes in the in- 
terpretation, and that the submillimeter transitions of HCO + 
and its isotope are tracing both outflowing material and the 
material that is participating in gravitational collapse near the 
cores towards both VLA 1623 and SM1N. 

3.2. Molecular Line Maps 
3.2.1. Outflows 

Outflows play an important role in the early phases of stellar 
evolution. In fact, bipolar outflows in Class sources, like 
VLA 1623, indicate that outflows arise in star formation at 
a very early stage. The collimated outflow observed toward 
VLA 1623 has bee n studied extensively since its discovery by 
lAndre etalJlll990l) . 

Figure [3] presents the newly obtained millimeter and sub- 
millimeter CO maps towards the p Oph A region. The left 
panel shows the FCRAO CO maps in the J=l— >0 transition. 
Curiously, the grey-scale 13 CO image shows a ridge of emis- 
sion that peaks about 3' west of the submillimeter continuum 
sources, SMI, SM1N, VLA 1623 and SM2. The 13 CO map 
also shows a well-defined edge of emission to the east close to 
the location of the nearby B star S 1 . It has been suggested that 
the SM1N, SMI and SM2 cores form an arc of filamentary 
structures centered on the B star SI, possibly indicating a link 
between these s ubmillimeter clumps and star SI (MAN98, 
iDi Francesco et all fe004l) '). The left panel of Figure [3] also 
shows the red-shifted and blue-shifted line wing emission in 
12 CO J=l — ^0 emission. Centered on VLA 1623, we can see 
the outflow, especially the more collimated blue-shifted lobe 
of the outflow. The red-shifted line wing emission appears 
more diffuse throughout the region without a clear signature 
representing bipolarity of the outflow. A second blue-shifted 
component is seen east and north of VLA 1623. 

The submillimeter CSO CO J=3^2 data have higher an- 
gular resolution and in this higher transition, excitation con- 
ditions are more favorable for picking up the strong line 
wing emission from outflows. The right panel of Figure [3] 
shows the integrated intensity maps from the CO J=3 — >2 line 
wing emission over a smaller region in p Oph A. Indeed, in 
this transition the highly collimated outflow of VLA 1623 is 
readily seen extending across the map from the northwest to 
the southeast. The superposition of the red and blue-shifted 
emission indicates the VLA 1623 flow is located almost in 
the plane of the sky. Also, the differing angles (s=s 20°) of 
the northwest red and blue outflow lobes is con sistent with 
the CO J=2 — ^1 observations from lAndre et alJ £l990) and 
And re^etafl fi99l . The north-west blue-shifted high veloc- 
ity lobe attributed to VLA 1623 may be a combination of an 
outflow from VLA 1623 more in line with the VLA 1623 red- 
shifted outflow lobe and a separate blue-shifted outflow orig- 
inating from some source north and east of VLA 1623. North 
of the VLA 1623 flow, we see two well-separated red-shifted 
and blue-shifted lobes. The midpoint of the line joining the 



peaks of red-shifted and blue-shifted emission of this flow 
north of VLA 1623 falls very close to SM1N. We suggest 
that we have detected a new bipolar outflow centered on the 
source SM1N. 

A second outflow was reported to have been discovered 
north of the VLA 1623 flow at approximately the position 
reported here, based on smal ler angular extent mil limeter in- 
terferometric measurements (Kam azaki et alJ | 200lL hereafter 
KSHK ), and JCMT CO J=3^2 observations jKamazaki et al. 
120031 hereafter KSHUK). However, both of these studies cov- 
ered a region only 1' centered on SM1N in right ascension. 
Our maps cover 2.5' around SM1N in RA. As a result of the 
smaller scale of their observations, KSHK and KSHUK claim 
the detection of a new northern flow (which they dubbed 'AN 
outflow") with a driving source that is a near-infrared star 
called GY30 (which is ~ 35" from SM1N). They also con- 
clude that the AN outflow has a blue-shifted lobe in the east, 
and a red-shifted lobe in the west, which is the same as the ori- 
entation of the VLA 1623 flow. A comparison of Figure 2 of 
KSHUK with Figure|3]of this paper indicates that their blue- 
shifted lobe coincides with the northern blue-shifted lobe seen 
in our work. But KSHUK do not observe the region where 
the eastern red-shifted lobe is seen. As a result, they have 
misidentified the source and orientation of the AN outflow. 
It appears from our present work that a new bipolar outflow 
is present north of VLA 1623, with a driving source, which 
is most likely SM1N, and whose orientation of blue and red- 
shifted lobes in the sky is opposite to that of VLA 1623. 

We have estimated the physical parameters of the newly dis- 
covered outflow towards SM1N using the CSO CO J=3^2 
data (see Table |3jl. The outflow from VLA 1623 is not fully 
mapped here, so we do not attempt to repeat the estimates 
of physical parameters o f the VLA 1623 outflow (see for e.g. 
Andre eTalJI 19901 H9 93). In order to estimate the masses of 
the outflowing gas, local thermodynamic equilibrium (LTE) 
conditions are assumed with an excitation temperature of 
T ex = 40 K JAndreet alJll990l) . The total 12 CO gas column 
density towards each lobe of the SM1N outflow is estimated 
by using the integrated brightness temperature under the line 
for velocities of -8 to 2 km s" 1 and 5 to 15 km s" 1 respectively 
for the blue and red-shifted lobes, and assuming that the 12 CO 
emission is optically thin. A fractional abundance of X( l2 CO) 
= 1 x 10~ 4 and a distance of 160 pc to p Oph are used. 

We then estimated dynamical timescale, momentum, ki- 
netic energy, and mechanical luminosity of the flow. The dy- 
namical timescale, tj was estimated from outflow length and 
velocity that was assumed to be constant. The differences 
between the systemic velocity and the central velocities of 
blueshifted (redshifted) maps were used as the outflow veloci- 
ties V of blueshifted (redshifted) components. The calculated 
values of the momentum P = MV, kinetic energy E=MV 2 /2, 
and mechanical luminosity L m = E/tj are summarized in Ta- 
ble [3] These quantities are not corrected for a possible incli- 
nation angle of the outflow. 

This new discovery of the SM1N outflow may indicate ei- 
ther this protostar is more evolved than previously believed 
or the onset of outflows occur almost simultaneously with the 
first indication of infall. 

3.2.2. Impact of Outflows 

With detected outflows in both sources, we can determine 
the extent to which the HCO + line profiles are affected by the 
molecular outflows mapped in CO. To do this, we made inte- 
grated intensity maps of the blue-shifted, line-core, and red- 
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TABLE 3 

Physical Parameters of SM1N Outflow 
(Based on CSO CO J=3->2) 



-200 

FIG. 3. — Outflows in the region. The left panel shows the 13 CO J=l— >0 line-core emission (in grey-scale) and the 12 CO J=l— +0 line-wing emission (in 
contours) observed with the FCRAO 14 m. The "CO emission is integrated over 1 to 6 km s , and the 12 CO over -3 to 2 km s (blue shifted - dashed contours), 
and 5 to 10 kms~'(red shifted - solid contours). The contour levels in K-kms -1 for the J=l— »0 data on the left are as follows: 15 to 35 by 1.5 ( 13 CO), 6 to 
25 by 1.5 ( 12 CO blue), 6 to 26 by 2 ( 12 CO red). The inset to the right shows the inner 190" X 150" region mapped in 12 CO J=3^2 emission. Blue shifted 
emission (-8 to 2 km s ; dashed lines) has contour levels from 8 to 30 by 2 Kkm s while red shifted emission (5 to 15 km s ) has contour levels from 8 to 30 
by 2 K-kms -1 . The star on the left panel shows the location of the nearby B star SI (see MAN98), while the diamonds in both panels show the location of the 
embedded submillimeter continuum objects. The square shows the location of the 2.111 fim near-infrared source GY30 (see KSHK). The X and Y axes of both 
plots show the angular offsets in arc-seconds from the position of VLA 1623. 

line-wing emission of HCO + confirms the detection of the 
outflow seen in CO J=3^2 emission. There are clear in- 
dications that the line-wing emission is elongated along the 
outflow direction. 

The HCO + line- wing integrated intensity about VLA 1623 
is more complex. There is a definite elongation of red and 
blue-shifted emission along the well known outflow driven by 
VLA 1623. Although there is a slight deviation at the edges of 
the map, the line-core, blue-shifted, and red-shifted emission 
are spatially coincident along the outflow. In addition to the 
outflow elongation, the red-shifted and blue-shifted emission 
has a small component perpendicular to the outflow. The red- 
shifted emission is predominately toward the north-north-east 
and the blue-shifted emission is stronger in the south-south- 
west direction. This emission perpendicular to the outflow 
direction may be tracing the overall rotation of the molecular 
cloud about VLA 1623. 



Parameter 


Blue 


Red 


Velocity (kms~') 


6 


5.5 


Length (xlO 4 AU) 


1.2 


1.1 


Mass(xl0~ 3 M ) 


5.5 


3.5 


Dynamical Timescale (x 10 3 yrs) 


9.2 


10 


Momentum (x 10 -2 Mokms -1 ) 


3.3 


1.9 


Kinetic Energy (x 10 42 ergs) 


2.0 


1.0 


Mechanical Luminosity (x 10~ 3 Lq) 


1.7 


0.8 



shifted emission of HCO + J=4^3. Figure|4|shows the HCO + 
J=4^3 integrated emission around VLA 1623 and SM1N, 
where the line core integrated intensity is shown in grayscale, 
blue-shifted emission in dashed contours, and the red-shifted 
emission is shown in solid contours. The line-core emission is 
expected to trace the denser, infalling material, while the line- 
wing maps may show contamination due to material swept-up 
by the outflows. 

Figure |4] shows a complex morphology of both the line- 
core and line-wing emissions. The line-core emission dis- 
plays two orthogonal components, one elongated along the 
direction of the outflows and one perpendicular. This effect 
is readily seen in the VLA 1623 region and to a lesser ex- 
tent in the region around SM1N. The overall structure of the 
line-core emission is consiste nt with the c ontinuum map of 
the cloud core region in Andre et al] i 19931) . Figure 1 . 

The line-wing emission around SM1N is collimated and 
consistent with the high velocity CO map in Figure [3] The 



3.2.3. Centroid Velocity Maps 

Given a velocity range, the centroid velocity is defined 
as that velocity at which the integrated intensity is equal 
on both sides. It has been shown that in the presence of 
complicated velocity fields, centroid velocity maps provide 
a better indication of the underlyin g velocity fields than inte- 
grated intensity maps jAdelson & Leundll9 88: Walk eTet alJ 
119941 iNaravanan & Walkerl ll998tlNaravanan et alJl2002t). In 
the model isovelocity maps of [Narayanan & Walker! (11991 
which considered a rotating infalling protostellar core, the ro- 
tational velocity field imposes a gradient of blue-shifted to 
red-shifted velocities, with the sense of the gradient being 
orthogonal to the rotational axis. When infall dominates in 
the central regions, the line profiles in the central region be- 
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FIG. 4.— CSO HCO + J=4^3 Integrated Intensity Maps toward VLA 1623 
and SM1N. Blue-shifted emission (-2 to 2 kms ) is shown in dashed con- 
tours, line core emission (2 to 6 km s~') is shown in grayscale, and red-shifted 
emission (6 to 10 km/s) is shown in solid contours. The contour levels are 0.1 
to 0.45 by 0.05 K-km/s for blue shifted emission, 2 to 6.5 by 0.5 K-km/s for 
line core emission, and 0.1 to 0.45 by 0.05 K-km/s for red-shifted emission. 

come blue asymmetric, and hence there is a preponderance 
of blue-shifted velocities in the central regions of the isove- 
locity maps, giving rise to the so-called "blue-bulge" signa- 
ture. Here, we make centroid velocity maps from our data 
to see if such an infall signature can be detected towards 
VLA 1623 andSMlN. 

We present the HCO + J=4^3 centroid velocity map of 
the VLA 1623 and SM1N region in Figure Blue-Shifted 
and red-shifted emission are shown with dashed and solid 
isovelocity contours respectively. It has been shown in ^3.2.21 
that outflows show more of an impact on the line-wings of 
HCO + emission. In an effort to lessen the impact of outflows 
on the centroid velocity map, the centroid velocities are com- 
puted only over the line core emission (2-6 km s" 1 ). 

The isovelocity contours around VLA 1623 are parallel to 
the direction of the outflow axis and show a gradient from the 
north-east to the south-west. The resulting gradient could be 
interpreted as a general rotation of the gas around VLA 1623 
with the south-western gas rotating toward the observer and 
the north-eastern gas rotating away from the observer. Some 
evidence for this direction of rotation is also seen in the HCO + 
integrated intensity map of Figure |4] The blue-shifted con- 
tours in VLA 1623 are even seen to extend north-east of the 
axis of rotation. The encroachment of the blue-shifted con- 
tours across the rotation axis displays the predicted "blue- 
bulge" infall signature. 

The region around SM1N shows a more complex morphol- 
ogy in the centroid velocity map. If the gradient in isovelocity 
contours towards SM1N is interpreted solely as that produced 
by rotation, one major feature is that the sense of rotation is 
counter to that of VLA 1623 with the blue-shifted gas north of 
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FIG. 5.— CSO HCO + J=4->3 centroid velocity map of VLA 1623 and 
SMIN. The centroid velocities only in the line core (2 to 6 km s~' ) are shown. 
The centroid velocities are compared with that of the optically thin isotope 
(~ 3.8 kms" 1 ). Blue-Shifted velocities ( < 3.8 kms" 1 ) are shown in dashed 
contours, and red-shifted velocities ( > 3.8 kms ) are shown in solid con- 
tours. Contour levels are 3.0 to 4.3 km s" 1 by 0.1 km s" 1 . 

the source and the red-shifted emission toward the south. An- 
other interesting feature is that unlike in VLA 1623 the isove- 
locity contours do not lie parallel to the outflow axis in SMIN. 
Toward the eastern edge of the map, the contours curve to the 
expected orientation if the contours were to align with the ro- 
tation axis. However, in the central regions, the gradient ap- 
pears to be in a north-north-west to south-south-east direction. 
Such a direction for the gradient in HCO + emission might re- 
sult if even its line-core in SMIN is tracing both outflow and 
infall, in which case the resultant velocity gradient would be 
oriented at an angle intermediate between the outflow direc- 
tion and rotational gradient direction. Although certainly less 
defined than the blue bulge found toward VLA 1623, SMIN's 
blue bulge is evident from the blue shifted shifted isovelocity 
contours extending from the north-west through the source to 
the south-east. 

4. DISCUSSION 

4.1. Infall Modelling 

To constrain the infall parameters in both sources, we mod- 
eled our sources using a semi-analytic infall model based 
on the Terebey, Sh u, and Casen (TSC) collapse treatment 
jTerebev et alll984l). For a more co mplete explanation of the 
mode l used, see iNaravananl lll997l) : iNaravanan et alJ d!998l 
2002). 

The TSC model requires the use of three parameters: the 
sound speed a, the rate of cloud rotation ft, and the infall time 
t (or infall radius r,„/ = at). After the model protostar has been 
generated, values for the temperature, density, and velocity 
along the prescribed lines of sight (LOSs) are generated for 



Kinematics in p Oph A 



(a) 



(b) 





20 

AR.A. (arcsec) 



-20 



20 

AR.A. (arcsec) 



-20 



FIG. 6. — Centroid velocity maps towards VLA 1623 for (a) observed (left), and (b) modelled (right) in the HCO + J=4— >3 transition. The contour intervals 
are 2.8 to 4.5 by 0.1 kms" 1 , and -.006 to 0.002 by 0.001 km s~> respectively. 



any given viewing angle. The LOSs are defined by two ge- 
ometric angles: a, the angle out of the plane of the sky, and 
ip, the angle in the plane of the sky. After the program gen- 
erates all LOSs, a separate radiative transfer code calculates 
molecular and continuum emission along all LOSs using the 
assumption of local thermal equilibrium (LTE). Many model 
runs were made varying the infall and geometric parameters 
each time. The synthetic observations were then convolved to 
the resolution of the telescopes used in the observations. The 
best fit models were determined by eye by constraining the 
model spectra to the observed HCO + and H 13 CO + spectra. 




FIG. 7.— TSC Model fits for VLA 1623. Observed HCO+ J=4^3 (his- 
tograms) and model HCO + J=4— >3 (solid lines) spectral mosaics are shown. 
The velocity extent displayed is to 8 km s , and the temperature scale is - 1 
to 5 K. 

The TSC model does not include outflows, while VLA 1623 
and SM1N both have vigorous yo ung ou tflow systems. How- 
ever, as has been discussed in ^3.2.21 while the line-wing 
emission of HCO + seems to trace the outflows, the line-core 
emission itself seems to be less affected by the presence of 
outflows. Indeed, for VLA 1623, the centroid velocity map 
shown in Figure|5]shows predominantly the impact of rotation 



and infall, both of which are well-treated by the TSC formal- 
ism. In our infall modeling, we only use the line-core HCO + 
to constrain the infall parameters. Even in this light, the infall 
parameters derived here should be used with some caution, as 
the models are for single isolated protostars, and not for clus- 
ters of protostars, and hence not rigorously applicable to the 
observations. 

The comparison between the observed spectra and the sim- 
ulated spectra of the best-fit model towards VLA 1623 is 
shown in Figure While the spectral shapes are reasonably 
well-fit (at least in the central regions), the best fit TSC model 
does not do a good job of matching the observed centroid ve- 
locity map observed in HCO + J=4-^3 (see Figure [6}. Given 
the caveats discussed above of applying the TSC model for 
this study, this latter result is not surprising. Similar model 
plots were made for SM1N, but are not shown here. 

Since there is some evidence that even the line-core emis- 
sion of HCO + inSMlN is being impacted by the observed 
outflow (see § 13.2.31 and Figure it was expected that it 
would be harder to fit the SM1N observations to the TSC 
model. The velocity gradient in Figure [5] seems to imply a 
more or less east-west rotational axis if the gradient is purely 
interpreted as due to that of rotation. In view of the pos- 
sible contamination from the outflow, in our modeling, we 
forced the rotational axis for SM1N to be parallel to that of 
VLA 1623. Because of this, the model centroid velocity map 
for SM1N cannot reproduce the direction of the gradients seen 
in the observed map of Figure|5] However, the general appear- 
ance and width of the line profiles are well modeled even in 
the region towards SM1N. 

When the effects of outflows are separable from the obser- 
vations of density-sensitive molecules such as HCO + and CS, 
the TSC model can be used reasonably effectively to derive 
infall parameters (see t he c ase of IRAS 16293 and SMM4 in 
iNaravanan et al.Nl 9981) and lNaravanan et alJ<2002ft '). The re- 
sults of the TSC models clearly indicate the difficulty in mod- 
eling complex regions that consist of both infalling and out- 
flowing motions. In fact, the inability to adequately account 
for the presence of outflows, severely handicaps the modelling 
of spectra and centroid velocity maps along the axes of the 
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TABLE 4 
Best Fit Model Parameters 



Source r inf a l inf Q a ip v turb 

(pc) (kms-') (yr) (s^ 1 ) (°) (°) (kms" 1 ) 

VLA 1623 0.01 0.33 3.0 x 10 4 3 x 10" 14 -60 0.5 

SM1N 0.006 0.58 l.OxlO 4 7xKT 14 -60 0.6 



outflows. Despite this shortcoming and the complexity of the 
region, our models are robust enough to deduce key properties 
of the regions around VLA 1623 and SM1N. The TSC mod- 
els confirm VLA 1623's outflow is within the plane of the sky 
to ±15° as initially proposed by Andre et al. ( 1990). Second, 
our modelling of VLA 1623 could not reproduce the obser- 
vations unless a temperature gradient was applied to the TSC 
code. The best fit temperature gradient decreases at an angle 
of 32° west of south. The nearby B3 star, SI, provides the 
best physical explanation for this temperature gradient. Our 
best model parameters are summarized in Table|4| 

4.2. Timescales 

One of the more difficult challenges in star formation is the 
determination of age. This difficulty arises from the lack of a 
definitive definition for the age of a protostar and inability to 
accurately measure key time dependent observables such as 
mass infall rate and velocity dispersion. One of the first es- 
timate s for the age of VLA 1623 was derived bv lAndre et alJ 
Jl990T) bv examining the dynamical timescale for VLA 1623's 
molecular outflow. Their conclusion, a dynamical age of 
-3000 yr, is similar to the age derived by AWB93 (<6000 
yr) using the mass infall rate. 

Using our submillimeter CO J=3 — >2 data, the dynamical 
age of VLA 1623 i s found to be —4000 yr, in agreement with 
lAndre et alJ d!990h . However, our modelling of VLA 1623 
suggests an older age for the time since infall by approxi- 
mately 10 times (see Table |4j. Given errors in the calcula- 
tions of both the dynamical and modelled age estimates, the 
similarity of the infall and dynamical timescales leads to the 
conclusion that large outflows develop almost from the onset 
of collapse. 

5. CONCLUSIONS 

We have observed two cores in the p Oph A molecular 
cloud to investigate the kinematic signatures of star forma- 
tion. We obtained isotopic and main line CO J=3— >2, J=l— *0, 



and HCO + J=4^3, J=l — >0 data of a known Class source, 
VLA 1623, and a core that was previously identified as a pre- 
protostellar core, SM1N. Using more spatially extended and 
sensitive CO maps than those of KSHK and KSHUK, we con- 
firm the detection of the second molecular outflow north of 
VLA 1623. This second outflow is most likely generated by 
the protostar SM1N, and not by the near-infrared star GY30 
as suggested by KSHK and KSHUK. SM1N appears to be in 
a more progressed evolutionary stage of star formation than 
previously thought. 

Millimeter observations toward the source positions show 
the expected blue asymmetry. However, the submillimeter 
transitions of HCO + J=4— >3 display a line profile signature of 
expanding rather than infalling gas. Most likely, in VLA 1623 
and SM1N, the higher density, centrally condensed gas kine- 
matics are greatly affected by the outflows emanating from 
both sources. Examination of the opacity profiles, however, 
provides support that, despite outflowing material, the core is 
still undergoing gravitational collapse. 

We performed centroid analysis of VLA 1623 and SM1N 
to confirm the existence of the "blue bulge" signature found 
in other star forming regions. The "blue bulge" found around 
VLA 1623 is very well aligned with the axis of rotation and 
is strong, confirming that VLA 1623 is a true protostar. The 
centroid velocity map near SM1N also exhibits "blue bulge" 
signature but it is not as well aligned about a rotation axis as 
seen in VLA 1623. 

The detection of infall and outflow motions towards SM1N 
suggests that both kinematic motions occur simultaneously 
from the earliest stages of star-formation. 
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